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Summary  of  Results 


A.  General  Accomplishments.  Progress  on  several  of  the  speci¬ 
fic  tasks  reported  below  deserved  special  emphasis.  These 
include  the  successful  execution  of  the  ARPA-CLIMAP  18,000 
YBP  climate  simulation  experiment  (discussed  as  Task  15)  ; 
the  results  obtained  during  the  conference  on  Transfer 
Functions  held  at  Madison,  Wisconsin,  April  3-5,  1974  (Task 
1) ;  the  assembly  of  pollen  maps  at  selected  times  during 
the  Holocene  (Task  9t) ;  the  progress  on  air-sea  interaction 
studies  represented  by  work  on  the  Gulf  Stream  off  Cape 
Hatteras  (Task  10a)  and  at  higher  latitudes  (Task  9v) ;  and 
the  development  of  a  new  uplift  model  for  reef  terrace 
sequences  (Task  8r) . 

In  addition,  there  are  two  ARPA-CLIMAP  related  activi¬ 
ties  which  are  time  consuming  and  do  not  appear  on  the  pre¬ 
scribed  list  of  tasks.  First,  as  a  member  of  the  GARP 
Panel  on  Climatic  Variation,  Imbrie  made  a  major  contribu¬ 
tion  to  those  portions  of  the  report  (Gates  iet  aJ. ,  1974) 
that  deal  with  paleoclimates .  Second,  as  a  member  of  the 
US-USSR  Joint  Committee  on  Cooperation  in  the  Field  of 
Environmental  Protection,  Imbrie  is  scheduled  to  spend  two 
weeks  in  the  Soviet  Union,  June  8-21,  1974.  He  and  J.  M. 
Mitchell,  Jr.  are  members  of  Sub-group  4  of  Working  Group 
8,  dealing  with  evidence  from  past  climates. 
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B.  Specific  progress  on  project  tasks.  On  pages  10-13  of 
the  Brown  ARPA  proposal  dated  July  18,  1972,  a  research 
plan  centered  on  sixteen  numbered  tasks  was  presented. 

These  tasks  are  listed  below.  Progress  made  towards  them 
is  discussed  in  brief  progress  statements  lettered  con¬ 
secutively  (a) .  (d*).  Those  tasks  numbered  without 

parentheses  will  be  carried  out  as  ARPA  research  and  will 
be  supported  exclusively  or  predominantly  by  ARPA  funds. 
Task  1:  TEST  ALTERNATE  TRANSFER  FUNCTION  METHODS 

(a)  A  Transfer  Function  Conference  was  held  at  Madi¬ 
son,  Wisconsin,  April  3-5,  1974,  jointly  funded  by 
ARPA  and  by  the  National  Science  Foundation's  IDOE 
Program.  The  organizing  committee  included  Imbrie 
and  Webb  of  Brown  University  (who  took  major  scienti¬ 
fic  responsibility) ,  Kutzbach  and  Wendland  (University 
of  Wisconsin) ,  Fritts  (University  of  Arizona) ,  and 
McIntyre  (Lamont-Doherty ) .  Approximately  30  scien¬ 
tists  attended  including  N.  Shackleton  and  H.  J.  B. 
Birks  (University  of  Cambridge)  and  T.  A.  Wijmstra 
(University  of  Amsterdam) .  A  total  of  15  institu¬ 
tions  and  organizations  was  represented  including: 
University  of  Wisconsin,  Brown  University,  Lamont- 
Doherty  Geological  Observatory  of  Columbia  University, 
Oregon  State  University,  University  of  Arizona,  Uni- 
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versity  of  Washington,  West  Georgia  College,  NOAA, 
University  of  Rhode  Island,  University  of  Minnesota, 
U.S.  Geological  Survey,  University  of  Miami,  U.S. 

Naval  Oceanographic  Office,  Rensselaer  Polytechnic 
Institute ,  and  Yale  University.  The  Conference  ob¬ 
jectives  --  to  establish  working  relationships  among 
workers  in  the  transfer  function  field;  to  evaluate 
the  relative  merits  of  competing  algebraic  methods; 
to  assess  the  reliability  of  climatic  estimates;  and 
to  determine  research  still  needed  —  were  accomplished 
To  guide  the  conference,  Imbrie  and  Webb  prepared  and 
distributed  a  working  paper  (Appendix  A) . 

(b)  A  permanent,  informal  organization  of  workers  in 
the  transfer-function  field  was  agreed  to  at  the  Madi¬ 
son  Conference.  Mr.  Douglas  Clark  of  Brown  University 
agreed  to  distribute  a  Transfer-Function  Newsletter 
that  will  serve  as  a  communications  net  among  active 
workers . 

(c)  A  working  group  of  the  Madison  Conference  headed 
by  John  Imbrie  dealt  with  the  "No  Analog"  problem: 

How  can  paleoc lima tic  estimates  be  validated  when 
the  ecosystem  structure  may  have  been  different  from 
today's?  For  the  marine  data  a  plan  was  formulated 
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using  the  18,000  YBP  samples  accumulated  as  part  of 
the  ARPA-CLIMAP  modelling  experiment.  One  test  in¬ 
volves  the  simultaneous  use  of  estimates  derived 
from  different  groups  of  organisms.  A  second  test 
involves  the  use  of  those  core  samples  in  which  it 
is  possible  to  obtain  an  independent  estimate  of  sea- 
surface  temperature  from  the  isotopic  technique  oi 
Shackleton  and  Opdyke  (1973) .  As  reported  by  Shack- 
leton,  only  in  high-deposition  rate  cores  is  it  pos¬ 
sible  to  obtain  reliable  paleotempera ture  estimates 
using  oxygen  isotopes.  Those  estimates  are  derived 
by  subtracting  the  values  obtained  for  benthonic 
fossils  from  the  values  obtained  for  planktonic  fos¬ 
sils.  High  deposition  rate  cores,  however,  consti¬ 
tute  only  about  20  percent  of  those  used  in  preparing 
the  18,000  YBP  map. 

(d)  A  second  working  group,  headed  by  T.  Webb, 
dealt  with  the  problems  of  establishing  and  applying 
transfer  functions  to  pollen  data  on  time  scales  of 
0-1000  years  and  0-100,000  years.  Because  transfer 
functions  were  first  applied  to  pollen  data  spanning 
10,000  years  (Webb  and  Bryson,  1972),  the  assumptions 
and  procedures  of  this  initial  work  must  be  modified 
when  working  with  data  on  other  time  scales. 


ARL'A- supported  work  at  the  University  of  Wis¬ 
consin  (directed  by  Drs .  J.  E.  Kutzbach  and  A.  Swain) 
has  produced  a  0-1000  year  pollen  record  from  annually 
laminated  lakes,  with  these  data,  methods  must  be 
found  to  filter  out  pollen  variation  from  fires,  soil 
differences,  and  human  disturbance,  in  order  for  the 
climatic  message  to  become  evident. 

Few  cores  contain  a  pollen  record  of  100,000 
years  or  more.  The  pollen  data  in  these  cores  are 
critical  in  coordinating  the  land  and  sea  records. 

Three  of  these  cores  are  located  in  the  Mediterranean 
area  where  human  disturbance  of  the  vegetation  dates 
back  over  thousands  of  years  and  where  few  samples  of 
the  modern  pollen  currently  exist.  The  group  agreed 
that  pollen  changes  in  the  fossil  record  are  much 
larger  than  the  changes  imposed  by  human  disturbance 
and  that  modern  samples  from  the  area  should  be  gather¬ 
ed.  T.  A.  Wijmstra,  who  has  worked  on  some  of  these 
cores,  agreed  to  collect  and  analyse  the  deeded  modern 

samples  from  initial  analyses  of  the  core  in  Macedonia 
(Wijmstra,  1969) . 

(e)  A  paper  prepared  by  T.  Webb  and  D.  Clark  was 
presented  to  show  empirically  the  differences  in  esti¬ 
mated  climatic  values  that  result  from  differences  in 


-  13  - 

the  numerical  techniques  used  to  derive  transfer 
functions.  Significant  differences  appear  in  the 
analysis  of  fossil  records,  e.g.  3°C  versus  6°C 
changes  in  temperature.  These  results  emphasize 
the  importance  of  a  thorough  examination  of  the  alge¬ 
braic,  statistical,  and  biological  features  of  each 
technique  in  order  to  identify  the  optimal  technique 
for  a  given  set  of  data. 


Task  2:  DEVELOP  OPTIMAL  TREND-SURFACE  TECHNIQUES 

(f)  No  work  has  been  done  since  that  reported  in 
Semi-Annual  Technical  Report  No.  1  in  which  we  des¬ 
cribed  a  contouring  program.  We  have  decided  that 
it  is  not  feasable  at  this  time  to  carry  on  with  the 
more  sophisticated  "optimal  trend-surface"  objective 
stated  in  our  initial  proposal.  This  objective  was 
to  develop  formal  procedures  for  drawing  contours 
with  limited  control  points  and  take  into  account 
physical  oceanographic  theory,  instead  of  using  the 
simple  interpolation  schemes  characteristic  of  auto¬ 
mated  contouring  programs.  Discussion  with  oceano¬ 
graphers  convinced  us  that  research  on  the  problem 
is  not  likely  to  yield  results  soon. 


Task  3:  COLLATE  MODERN  POLLEN  DATA  FOR  EASTERN  U.S. 


(g)  Data  on  115  samples  of  modern  pollen  were  added 
to  the  800  samples  compiled  in  the  first  year.  The 
new  data  included  samples  from  northwestern  Nebraska, 
northern  Montana,  central  Florida,  and  central  Ontar¬ 
io.  These  data  were  derived  from  published  tables 
and  from  unpublished  tables  made  available  by  Dr.  J . 
H.  McAndrews  (Royal  Ontario  Museum,  Toronto) . 

(h)  Maps  were  prepared  showing  the  distribution  of 
each  of  the  major  pollen  types  in  central  and  eastern 
North  America  (Figures  1-7) .  These  maps  were  pre¬ 
sented  and  discussed  in  two  papers  at  the  INQUA  meet¬ 
ings  in  New  Zealand  in  December,  1973  (Appendix  B  and 
C) .  Manuscripts  describing  this  work  are  now  in  pre¬ 
paration  (Davis  and  Webb,  in  prep,  and  Webb  and 
McAndrews,  in  prep.) .  These  manuscripts  and  maps 
will  inform  scholars  of  the  locations  of  modern  sam¬ 
ples  and  will  illustrate  those  areas  most  in  need 

of  additional  sampling. 

(i)  As  part  of  the  analysis  and  presentation  of  pat¬ 
terns  in  the  distribution  of  pollen  data,  multivari¬ 
ate  statistical  techniques  are  important  (see  Webb, 
1974) .  These  techniques  reveal  which  patterns  are 
most  prominent  within  the  data  and  summarize  the  im- 


portant  information.  Because  several  possible  tech¬ 
niques  can  be  used,  each  utilizing  different  assump¬ 
tions  about  the  data,  work  was  undertaken  to  compare 
the  results  of  three  of  these  methods  —  e.g.,  prin¬ 
cipal  components  analysis,  canonical  variate  analysis 
(closely  related  to  multiple  discriminant  analysis) , 
and  cluster  analysis  —  on  a  subset  of  the  large  set 
of  modern  samples.  The  result  of  this  work  shows 
that  the  same  basic  structure  in  the  data  appears  no 
matter  what  technique  is  used.  Slight  differences 
in  the  details  of  each  technique,  however,  are  impor¬ 
tant  in  revealing  different  aspects  of  the  data.  A 
manuscript  is  now  in  the  final  stages  of  preparation 
describing  this  work  (Birks,  _et  _al .  ,  in  prep.). 

(j)  In  cooperative  work  at  the  University  of  Wiscon¬ 
sin,  Webb  provided  advice  to  Mr.  G.  Peterson  for  the 
compilation  of  modern  pollen  samples  from  the  Soviet 
Union. 

(k)  Close  working  ties  were  established  with  H.  J. 

B.  Birks  (Cambridge  Univ.,  England)  and  J.  H.  McAndrews 
(Royal  Ontario  Museum,  Toronto,  Canada)  for  data  ex¬ 
changes  and  for  sharing  methods  for  statistical  ana¬ 
lysis  of  data  (Birks  ejt  al_. ,  in  prep.)  . 
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Task  4:  IMPROVE  ABSOLUTE  CHRONOLOGY  OF  THE  30,000-130,000 
YBP  INTERVAL 

(1)  There  is  no  circa  40,000  YBP  terrace  on  Barba¬ 
dos.  A  terrace  approximately  two  meters  above  sea 
level  and  of  limited  geographic  distribution  was  pre¬ 
viously  suspected  to  be  circa  40,000  YBP.  Detailed 
mapping,  sampling,  and  carbon- 14  dating  now  clearly 
indicate  this  feature  to  be  less  than  4,000  years  old, 
essentially  Recent.  We  are  lead  to  the  firm  conclu¬ 
sion  that  any  high  stand  of  the  sea  which  occurred 
in  the  time  range  10,000  to  60,000  BP  did  not  reach 
high  enough  to  be  recorded  above  sea  level  on  Barbados 

(m)  Additional  dates  are  being  obtained  on  the 
"60,000  YBP"  Barbados  terrace.  The  60,000  YBP  inter - 
stadial  is  recorded  in  temperature  estimates  from 
many  North  Atlantic  cores  and  appears  to  be  the  last 
major  warm  interval  before  the  onset  of  full  glacial 
conditions.  It  is  therefore  important  that  we  get  as 
accurate  an  estimate  as  possible  of  the  timing  and 
elevation  of  this  high  stand.  Additional  Barbados 
samples  are  being  analyzed  and  will  hopefully  accom¬ 
plish  this  purpose. 

Task  5 :  EXTRACT  POLLEN  DATA  FROM  THE  LITERATURE  AND  INTER 
PRET  IT  IN  TERMS  OF  CLIMATIC  TIME  SERIES 
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(n)  The  data  obtained  from  50  long  pollen  cores 
with  radiocarbon  dates  from  0  to  15,000  B.P.  have 
been  key  punched  and  stored  on  punch  cards.  in 
addition,  data  from  20  other  long-cores  have  been 
added  to  this  data  set.  All  the  time-stratigraphic 
information  (mainly  radiocarbon-dates)  have  been 
recorded  and  is  being  used  in  selecting  pollen  data 
at  1,000  year  time-intervals  during  the  Holocene 
for  presentation  on  maps  (see  Task  9) . 

Task  6:  OBTAIN  AND  INTERPRET  NEW  POLLEN  DATA 

(o)  In  light  of  these  initial  maps,  new  samples 
have  been  collected  in  central  Nebraska  through 
cooperation  with  Dr.  L.  V.  Benson  (Dept,  of  Geology, 
Univ.  of  Nebrasks).  These  samples  provide  data  from 
previously  unsampled  areas.  Data  from  Kentucky, 
Illinois,  and  Missouri  have  also  been  collected, 
processed,  and  added  to  the  data  set  through  coopera¬ 
tion  with  palynologis ts  at  the  University  of  Wisconsin. 

Task  (7):  OBTAIN  SEA-SURFACE  PALEOTEMPERATURE  DATA  FROM 
HIGH-DEPOSITION-RATE  MARINE  CORES 

(p)  Work  is  proceeding  on  high-deposition  rate  cores 

in  the  Curiaco  Trench  (off  Venezuela) .  Although  we 

have  identified  a  180  and  380  year  periodicity  in 

the  faunal  data,  reliable  paleoclimatic  estimates  are 


. .  .  . 


Task 


not  yet  possible  because  of  the  restricted  fauna. 

We  are  currently  preparing  a  data  set  from  which  we 
can  derive  a  more  reliable  transfer  function  for  this 
unusual  fauna. 

Additional  high  deposition  rate  cores  have  been 
identified  from  the  Santa  Barbara  basin  off  Califor¬ 
nia.  Those  cores  have  varves  that  can  be  correlated 
with  varves  in  dated  cores  back  to  about  1100  A.D. 

Work  is  in  progress  on  extending  these  correlations, 
and  analyzing  the  fauna  in  these  cores. 

8:  DETERMINATION  OF  SEA-LEVEL  FLUCTUATIONS  IN  THE 

RANGE  75,000-130,000  YBP 

(q)  Progress  on  this  task  in  general  is  reported  by 
Matthews  (1974)  (abstract  attached  as  Appendix  D) 
and  Bloom  et.  al .  (1974)  (abstract  and  key  figures 

attached  as  Appendix  E)  .  Barbados  and  New  Guinea 
data  are  compatible  to  a  first  approximation  and  in¬ 
dicate  numerous  high  stands  of  the  sea  separated  by 
low  stands  of  considerable  magnitude  (20  to  70  meter 
variations  in  sea  level  within  5,000  to  10,000  years). 

(r)  A  new  uplift  model  has  been  developed  (Appendix 
F)  that  affords  a  more  precise  evaluation  of  tectonic 
effects  on  terrace  sequences  over  short  intervals  of 
RKM  is  presently  collaborating  with  John  Chap- 
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pell  of  Australian  National  University  concerning 
application  of  the  model  to  New  Guinea  data.  This 
study  should  produce  more  precise  estimates  of  the 
history  of  sea-level  dynamics. 

(s)  The  125,000  YBP  high  stand  lasted  no  longer 

than  5,000  years,  was  pr  ,ded  by  a  25  meter  sea 

level  rise,  and  was  followed  by  sea  level  lowering 

at  a  rate  of  5  to  10  meters  pe  thousand  years. 

Mangion  and  Matthews  (1974)  abstract  is  attached  as 
Appendix  G. 


^POLLEN  T“E  SERIES  F°R  N°RTH 

(t)  Forty- four  maps  have  been  produced  to  show  the 
distribution  of  major  pollen  types  (spruce,  pine, 
oak,  and  herbs)  at  1000  year  intervals  during  the 
Holocene.  In  addition,  difference  raps  ha'  e  been 
constructed  to  show  changes  in  pollen  distributions 
from  one  time  interval  to  the  next.  These  maps 
(Fig.  8)  show  clearly  the  changes  in  the  pollen  (and 
by  inference  vegetation  and  climate)  as  the  ice  re¬ 
treated  northward.  A  manuscript  (Bernabo,  Webb,  and 
McAndrews,  in  prep.)  is  now  being  prepared  to  pre¬ 
sent  these  maps.  This  manuscript  will  illustrate  to 
paly.iologists  where  more  samples  are  needed,  espe- 
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cially  if  synoptic  maps  of  the  pollen  data  are  to 
be  constructed.  Although  mapping  the  data  is  criti¬ 
cal  for  paleoclimatological  studies,  few  Quaternary 
palynologists  have  attempted  to  map  their  data  (see 
McAndrews  and  Power,  1973;  Webb,  1974) . 

(u)  Work  continued  on  finding  or  writing  computer 
programs  for  compiling,  editing  and  updating  pollen 
data.  Work  was  begun  to  tie  these  file-handling 
programs  in  with  the  programs  for  finding  transfer 
functions  and  for  plotting  or  mapping  the  output. 

M  A  project  was  started  with  Dr.  W.  Ruddiman  (U.S. 
Naval  Oceanographic  Office)  to  study  the  teleconnec¬ 
tions  between  climatic  patterns  over  eastern  North 
America  and  the  central  North  Atlantic.  This  pro¬ 
ject  combines  the  paleoclimatic  information  of  pollen 
data  with  the  paleo-oceanographic  data  of  planktonic 
foraminifera . 


Task  10: 


OBTAIN  NEW  CLIMATIC  TIME  SERIES 
RATE  MARINE  CORES 


IN  HIGH-DEPOSITION 


(w)  Two  high  deposition  rate  cores  have  been  identi¬ 
fied  off  Cape  Hatteras,  North  Carolina.  Sedimenta¬ 
tion  rates  between  radioes:: bon-dated  postglacial  and 
late  glacial  horizons  reach  27cm/1000  years;  the  full 
glacial  sedimentation  rate  probably  exceeded  50cm/1000 
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years.  By  examining  pollen  and  foraminif era  from 
identical  core  samples,  oceanographic  changes,  in 
particular  Gulf  Stream  migrations,  can  be  directly 
correlated  with  changes  in  terrestrial  vegecation. 
Initial  results  indicate  that  there  is  a  slight  time 
lag  (500  to  2000  years)  between  changes  on  land  and. 
sea,  with  the  sea  leading  the  land.  A  manuscript 
describing  this  work  is  being  prepared  (Balsam  and 
Florer,  in  preparation) . 

(x)  We  have  identified  high-deposition-rate  cores 
in  the  eastern  Mediterranean  that  are  amenable  to 
transfer  function  analysis.  Data  presented  in  papers 
by  Todd  (1958)  and  by  Parker  (1958)  have  been  analyzed 
in  our  laboratory  by  transfer  function  techniques. 

The  result  is  an  equation  which  satisfactorily  esti¬ 
mates  pal eotempera tures  and  salinity  in  Mediterranean 
cores.  Sufficient  C-14  and  stratigraphic  data  are 
available  to  prove  that  at  least  one  core  off  the 
Nile  delta  (Core  184)  has  a  deposition  rate  of  about 
40cm/1000  years,  and  penetrates  the  entire  Holocene. 
This  opportunity  will  be  exploited  during  the  coming 
year . 

Task  (11) :  CORRELATE  CLIMATIC  RECORD  OF  THE  ICE  CORES 


(y)  A  short  paper  (Andrews  e_t  _al_.  ,  1974,  in  press) 

represents  the  application  of  correlations  made  by 
Sancetta  et  aT . ,  (1973)  be  tween  North  Atlantic  core 

V23-82  and  the  Camp  Century  ice  core  record.  Andrews 
et  al . ,  conclude  that  as  suggested  by  Imbrie  (1972) 
snowfall  in  the  Baffin  Island  region  was  lower  than 
normal  during  the  18,000  YBP  maximum-glaciation  in¬ 
terval  . 

Task  (12) :  COOPERATE  WITH  GLACIAL  GEOLOGISTS  IN  OBTAIN¬ 
ING  ACCURATE  ICE  MARGIN  POSITIONS  FOR  SELECTED 
TIMES  DURING  THE  PAST  130,000  YEARS 

(z)  A  critical  element  in  the  ARPA-CLIMAP  18,000  YBP 
modelling  experiment  was  the  availability  of  accurate 
data  on  the  margins  and  elevations  of  ice-sheets. 
These  data  were  assembled  on  a  global  basis  by  Dr. 
George  Denton  of  the  University  of  Maine  and  Dr. 

Bjorn  Anderson  of  Norway.  This  compilation  was  made 
expressly  for  the  ARPA-CLIMAP  experiment.  A  confer¬ 
ence  held  at  Brown  during  December,  1973,  made  the 
final  integration  between  ice-margin  and  sea-surface 
temperature  data  used  for  the  experiment. 

Task  (13) :  IDENTIFY  SIGNIFICANT  FREQUENCIES  IN  CLIMATIC 
TIME  SERIES 

(a1)  Research  is  in  progress  on  long  Pacific  deep- 
sea  core  records,  in  cooperation  with  N.  Shackleton, 


of  Cambridge.  The  initial  result  w as  a  spectrum 
for  Pacific  core  V28-238,  which  showed  a  significant 
100,000-year  periodicity  and  insignificant  spectral 
peaks  in  the  vicinity  of  40,000  and  20,000  years 
(Fig.  2  of  Semi-Annual  Technical  Report  No.  1  was 
reproduced  in  Appendix  A  of  the  Report  to  the  GARP 
Panel  on  Climatic  Change.)  Current  work  will  result 
in  a  spectrum  from  a  second  long  Pacific  core  (V28- 
239) .  If  the  same  frequencies  appear,  we  will  pub¬ 
lish  the  results. 

Task  (14) :  STUDY  REGIONAL  VARIATION  IN  CLIMATIC  TIME 
SERIES 

(b1)  This  work  has  led  to  one  important  discovery: 
that  significant  phase  differences  occur  in  the  re¬ 
sponse  times  of  various  parts  of  the  air-sea-ice 
system  during  the  glacial-interglacial  cycle.  In 
particular,  many  parts  of  the  tropical  and  sub¬ 
tropical  sea-surface  achieve  their  minimum  tempera¬ 
tures  several  thousand  years  before  18,D00  YBP,  when 
both  the  Northern  Hemisphere  ice-sheets  and  the 
oceanic  polar  front  in  the  North  Atlantic  achieve 
maximum  southward  extent.  This  fact,  reported  by 
Hays  (1974),  was  interpreted  by  Imbrie  (1974)  as  con¬ 
tradicting  the  standard  Milankovitch  theory  of  the 
ice  ages.  For,  if  the  ice-age  is  triggered  by  a 
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high-latitude  snow-line  reaction,  then  the  low- 
latitude  ocean  should  respond  in  phase  with  or  lag 
behind  the  ice-sheets.  This  line  of  argument  has  led 
Imbrie  to  a  re-casting  ol  the  theory  of  orbital  con¬ 
trol  of  climate,  stressing  low-latitude  effects.  Work 
along  this  line  forms  a  significant  part  of  our  plan 
for  research  during  the  next  two  years . 

Task  15:  COOPERATE  WITH  ARPA  SCIENTISTS  IN  THE  DESIGN  AND 
EXECUTION  OF  EXPERIMENTS  TO  TEST  NUMERICAL  CLI¬ 
MATE  MODELS  AGAINST  SYNOPTIC  PALEOCLIMATIC  DATA 

(c  )  Work  on  this  task  has  been  a  major  pre-occupation 

of  the  Brown-ARPA  group  during  the  past  year.  These 

of for ts  came  to  fruition  on  April  8,  1974,  when  an 

eight-paper  session  at  the  Meteorology  Section  of 

the  AGU  Meeting  in  Washington,  D.C.,  was  devoted  to 

a  description  of  the  design  and  execution  of  the 

ARPA-CLIMAP  August  18,000  YBP  simulation  experiment. 

This  set  of  abstracts  is  reproduced  as  Appendix  H. 

Gates  reported  a  sea-level  pressure  pattern  calculated 

for  August  18,000  YBP.  Preliminary  analysis  by  Gates 

indicates  that  an  intuitively  reasonable  result  was 

obtained .  Current  ARPA— CLIMAP  activity  is  focussed 

on  the  use  of  independent  data  to  check  the  accuracy 

of  the  model.  The  bulk  of  this  test  data  represent 

pollen  spectra  from  around  the  world,  compiled  and 
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ancxlyzed  mainly  at  Brown.  Conferences  are  planned 
during  late  June  and  September  in  which  the  testing 
and  evaluation  of  the  RAND  climatic  model  will  be 
carried  out. 


Task  16:  IDENTIFY  SYNOPTIC  PATTERNS  OF  THE  PAST  20,000 

YEARS  WHICH  ARE  ANALOGOUS  TO  EXTREMES  DURING  THE 
PAST  130,000  YEARS 

(d1)  Work  on  this  task  will  begin  when  additional 
maps  depicting  the  past  20,000  years  and  the  past 
130,000  years  become  available.  We  were  premature 
in  proposing  to  accomplish  this  task,  when  we  wrote 
the  initial  proposal  in  1972. 
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IV.  TECHNICAL  REPORT  SUMMARY 

A.  Purpose  of  the  research.  Climate  is  always  changing. 

The  practical  as  well  as  the  purely  scientific  value  of  understand¬ 
ing  the  processes  which  bring  about  climatic  change  is  self-evident. 
Only  by  understanding  these  processes  can  we  comprehend  past  and 
predict  future  climates.  To  date  this  goal  has  not  been  achieved. 
What  sort  of  research  is  needed?  In  part,  the  answer  will  come  from 
an  improved  understanding  of  the  mechanisms  by  which  the  global  air- 
sea-ice  system  yield  a  climate  in  equilibrium  with  today's  boundary 
conditions.  Significant  progress  in  this  direction  has  been  made 
by  several  numerical  models  of  climate.  But  this  research  cannot 
solve  the  puzzle  of  climatic  change  by  itself.  We  must,  in  addition, 
determine  what  changes  in  boundary  conditions  (if  any)  force  clima¬ 
tic  change  —  and  understand  the  forcing  functions  themselves. 

These  objectives  can  only  be  achieved  by  studying  the  workings  of 
the  global  climate  machine  over  a  time  span  adequate  to  record  a 
representative  range  of  conditions  in  Nature's  own  laboratory.  The 
purpose  of  ARPA  paleoclimate  research  at  Brown  is  to  document  and 
understand  climatic  changes  on  time  scales  ranging  from  thousands 
to  hundreds  of  thousands  of  years . 

B.  Methodology.  Part  of  our  research  effort  aims  at  two 
basic  observational  problems  in  paleoclimatography :  1)  obtaining 

accurate,  quantitative,  paleoclimatic  estimates;  and  2)  constructing 
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an  absolute  time  scale.  The  first  problem  is  approached  by  eval¬ 
uating  and  testing  a  number  of  multivariate  transfer  function  methods 
applied  to  fossil  pollen  records  (obtained  from  bogs  and  lakes) ,  and 
to  fossil  plankton  (obtained  in  deep-sea  cores) .  The  second  problem 
is  approached  by  radiometric  dating  of  the  sea-level  record  contained 
in  the  uplifted  Pleistocene  coral  reefs  of  Barbados,  by  radiometric 
and  dates  on  marine  and  terrestrial  cores,  and  by  using  oxygen 

isotope  data  to  correlate  between  cores. 

Other  research  is  directed  towards  obtaining  synoptic  maps  of 
past  climates  that  illustrate  the  geographic  patterns  of  climatic 
change.  Pollen-based  maps  of  terrestrial  climates  are  largely  the 
product  of  ARPA  research.  Maps  of  past  marine  climates  are  obtained 
by  ARPA  investigators  working  in  cooperation  with  an  NSF-funded 
CLIMAP  project.  The  interpretation  of  these  data  is  discussed  below 
(Global  Atmospheric  Modelling  Experiment) . 

C .  Technical  Results 

1.  The  ARPA-CLIMAP  global  atmospheric  modelling  experiment 
has  been  run  (Task  15) .  ARPA  and  CLIMAP  scientists  used 
geological  data  to  establish  the  boundary  conditions  for 
a  full  glacial  period  (18,000  YBP) .  The  boundary  condi¬ 
tions  they  specified  were:  sea-surface  temperatures, 
sea-levels,  ice-margins,  and  surface  albedo.  This  body 
of  data  is  the  largest  and  most  diverse  body  of  paleo- 
climatic  data  ever  assembled  for  a  single  time  period. 


It  will  probably  be  refined  and  may  be  used  to  test 
second  generation  numerical  climatic  models.  Results 
of  the  experiment  will  be  validated  using  both  pollen 
data  gathered  by  Brown  personnel  and  marine  data  not 
used  in  running  the  experiment. 

The  initial  results  of  this  experiment,,  a  sea 
level  atmospheric  pressure  map  for  full  glacial  condi¬ 
tions  (18,000  YBP) ,  were  presented  at  the  A.G.U.  meet¬ 
ing  April  8,  1974.  Analysis  of  the  results  and  the 
validation  procedure  are  still  in  progress,  and  it  will 
probably  be  several  months  before  we  understand  some 
of  the  ramifications  of  the  experiment.  Howev.r,  our 
results  in  general  are  positive,  and  we  are  now  one 
step  closer  to  understanding  climatic  change  and  pre¬ 
dicting  future  climates. 

Participants  at  the  transfer  function  conference  made 
significant  progress  toward  understanding  how  various 
mathematical  techniques  modify  results  (Task  1) .  As 
a  result,  improved  transfer  functions  should  produce 
more  accurate  paleoclimatic  estimates.  Computer  pro¬ 
grams  for  establishing  the  transfer  functions  will 
also  be  more  widely  shared. 

The  development  of  a  new  uplift  model  for  reef  terrace 
sequences  will  enable  us  to  obtain  much  more  accurate 
figures  on  the  rate  of  sea  level  fluctuations  (Task  8) 
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Data  on  the  magnitude  and  rate  of  sea  level  fluctua¬ 
tions  (and  its  inverse,  change  in  glacial  ice  volume) 
are  critical  for  testing  theories  of  climatic  change. 

4.  The  geographic  range  of  our  pollen  data  set  has  been 
extended  and  arrangements  made  for  obtaining  data  from 
most  of  the  Northern  Hemisphere  (Tasks  3  and  5)  .  Sig¬ 
nificant  progress  has  been  made  filling  in  gaps  in  the 
modern  data  as  described  in  Task  6.  Maps  are  now 
available  showing  the  distribution  of  all  the  import¬ 
ant  pollen  types  in  eastern  North  America  (Task  3) . 

5.  Synoptic  maps  of  the  past  distribution  of  certain  pol¬ 
len  types  have  been  prepared  from  eastern  North  America 
for  selected  intervals  of  the  Holocene  and  Late  Glacial. 
These  maps  contain  the  first  clear  picture  of  how  the 
vegetation  in  a  region  responds  to  a  major  climatic 
change  and  what  the  rate  of  this  change  is.  The  maps 
illustrate  the  distribution  of  the  raw  data  from  which 
paleoclimatic  maps  will  be  produced  (Task  9) . 

6.  Initial  work  has  been  carried  out  toward  developing  a 
system  of  computer  programs  capable  of  handling,  trans¬ 
forming,  and  mapping  the  large  sets  of  geographically 
distributed  fossil  data  used  in  paleoclimatic  recon¬ 


struction  (Task  9)  . 


30 


Figure  1 


Distribution  of  pollen  percentages  xn 
sediment  samples  of  modern  pollen  for 
oyperaceae  (sedge) ,  Salxx  (willow) , 

Be  tula  (birch)  .  and  Almjs  (alder) 
pollen.  Percentages  are  based  on  a 
sum  of  total  pollen  counted  excluding 
aquatic  pollen  types.  Isopleths  are 
of  equal  pollen  percentage. 
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Figure  2  Distribution  of  pollen  percentages  in 
sediment  samples  of  modern  pollen  for 
Picea  (spruce) ,  Pinus  (pine) ,  Abies 
(fir) ,  and  Juniperus-Thu ja  (juniper- 
arborvitae)  pollen.  Percentages  are 
based  on  a  sum  of  total  pollen  counted 
excluding  aquatic  pollen  types.  Iso- 
pleths  are  contours  of  equal  pollen 
percentage. 
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Distribution  of  the  pollen  percentages 
in  sediment  samples  of  the  modern  pollen 
for  Tsuga  (hemlock) ,  Fagus  (beech) ,  Acer 
(maple) ,  and  Ulmus  (elm)  pollen.  Per¬ 
centages  are  based  on  a  sum  of  total  pol¬ 
len  counted  excluding  aquatic  pollen  types, 
Isopleths  are  contours  of  equal  pollen  per¬ 


centage  . 
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Figure  4 


Distribution  of  the  pollen  percentages 
in  sediment  samples  of  the  modern  pollen 
for  Quercus  (oak) ,  Carya  (hickory) ,  Ostrya- 
Carpinus  (hop-hornbeam)  ,  and  Fraxinus  (ash) 
pollen.  Percentages  are  based  on  a  sum  of 
total  pollen  counted  excluding  aquatic 
pollen  types.  Isopleths  are  contours  of 
equal  pollen  percentage. 
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Figure  5  Distribution  of  the  pollen  percentages 

in  sediment  samples  of  the  modern  pollen 
for  Populus  (poplar) ,  Plantago  and  Rumex 
(plantain  and  sorrel) ,  Ambrosia  (ragweed) 
and  Gramineae  (grass)  pollen.  Percentages 
are  based  on  a  sum  of  total  pollen  counted 
excluding  aquatic  pollen  types.  Isopleths 
are  contours  of  equal  pollen  percentage. 
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Figure  6  Distribution  of  the  pollen  percentages  in 
sediment  samples  of  the  modern  pollen  for 
Compos itae  (other-composite) ,  Chenopodiineae 
(pigweed) ,  Artemisia  (sage) ,  and  prairie  herb 
(grass,  pigweed,  other-composite,  and  sage) 
pollen.  Percentages  are  based  on  a  sum 
total  pollen  counted  excluding  aquatic  pol¬ 
len  types.  Isopleths  are  contours  of  equal 
pollen  percentage. 
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Figure  7  Distribution  of  the  scores  of  the  first 
(P.C.  1),  second  (P.C.  2)  third  (P.C.  3) 
and  fourth  (P.C.  4)  principal  components. 


44 


Figure  8  The  distribution  of  pine  pollen  at 

selected  times  during  the  deglaciation 
of  Eastern  North  America  (Bernabo  _et  al. , 
1974) .  Contours  are  lines  of  pollen 
frequency,  expressed  as  a  percent  of 
total  pollen.  Control  points  representing 
radiocarbon  dated  cores  are  indicated  by 
the  dots.  The  approximate  margi.  j  of 
the  Laurentide  ice-sheet  are  indicated 
by  the  stippled  pattern  (after  Bryson 
et  al. ,  1569) . 
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WORKING  PAPER  FOR  THE  CONFERENCE  ON  TRANSFER  FUNCTIONS  SPONSORED  BY 
CLTMAr  AND  ARPA  AND  HELD  AT  THE  UNITES  Il'Y  OF  WISCONSIN,  April  1-5,  1^74 

John  Imbrie  and  Thompson  Webb  III 

I .  PREFACE . 

Included  in  this  working  paper  are  definitions  of  key  concepts 
and  a  list  of  problems.  The  definitions  are  subject  to  change  and 
the  list  is  incomplete.  The  aim  of  the  paper  is  to  improve  communi¬ 
cation;  to  provide  an  initial  frame  of  reference  within  which  to 
carry  on  debate;  and  to  indicate  the  scope  of  the  conference.  We 
hope  that  it  will  aid  the  conference  in  achieving  one  of  its  goals: 
the  design  of  experiments  to  answer  basic  questions  in  the  field  of 
transfer  functions. 

It  would  be  a  simple  matter  to  broaden  the  conference,  not  only 
in  the  field  of  biologically  based  transfer  functions,  but  also  in 
the  field  of  physically  based  transfer  functions.  Left  out  of  con¬ 
sideration  are  taxonomically  based  work  on  benthic  organisms,  amino- 
acid  techniques,  and  various  approaches  in  the  fields  of  glacial  and 
peri-glacial  geology,  to  name  only  a  few.  Both  in  writing  this  paper 
and  organizing  the  conference,  we  have  stressed  a  group  of  extensively 
tested  techniques  which  are  based  on  an  algebraically  homogeneous  set 
of  models  and  which  provide  continuous  time-series  and  synoptic  maps 
of  late  Pleistocene  climates.  Emphasis  is  placed  first  on  techniques 
which  monitor  biological  responses  to  climatic  change  (e.g.  Imbrie 
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and  Kipp,  1971;  Webb  and  Bryson,  1972,  Fritts  et  al . ,  1971);  and 
second  on  the  model  developed  by  Shackleton  and  Opdyke  (1973)  which 
uses  oxygen-isotope  data  from  benthic  and  planktonic  fossils  to 
Provide  physically  based  estimates  of  changes  in  surface-water  tem¬ 
peratures. 

The  mapping-potential  of  these  Lechniques  is  particularly 
important,  as  noted  in  the  Report  of  the  May  1973  meeting  of  SCOR 
Working  Group  40.  The  initiative  for  the  Wisconsin  Conference  stems 
from  this  Report,  in  which  Chairman  Van  Andel's  summary  states  that 
there  has  been  a  "rapid  development  in  recent  years  of  mathematical 
and  statistical  techniques  that  have  for  the  first  time  made  it 
possible  to  quantify  environmental  parameters  derived  from  the 
oceanic  geologic  record  in  a  form  suitable  for  input  to  dynamical 
models  of  the  circulation.  By  providing  quantitative  environmental 
data  for  past  conditions,  it  is  not  only  becoming  possible  to  model 
the  dynamics  of  past  circulations,  but  also  to  test  existing  models 
the  present  circulation  with  regard  ♦■o  their  power  of  explanation." 
The  Report  goes  on  to  recommend  the  organization  of  "a  small  techni¬ 
cal  workshop  on  mathematical  and  statistical  models  that  permit 
quantification  of  environmental  data  from  the  geologic  record... 
and  that  the  purpose  of  this  workshop  to  review  the  various  models, 
improve  their  compatibility,  assess  their  strengths  and  weaknesses, 
evaluate  the  need  for  further  methodological  work,  and  produce  a 
report  that  promotes  wider  use  and  understanding  of  these  techniques." 


II .  DEFINITIONS . 
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Let  the  matrix  X  be  a  defined  set  of  response  properties,  biotic 
or  isotopic,  measured  over  a  defined  realm  of  time  and  space.  Let  C 
be  a  set  of  physical  parameters  of  clrmate,  either  marine  or  atmos¬ 
pheric,  measured  over  the  same  time-space  realm  and  assumed  to  be 
causally  related  to  X.  Let  D  be  a  set  of  other  physical  parameters 
of  the  system  which  are  independent  of  C  and  which,  together  with  C, 
completely  control  the  response  X.  D  includes  such  factors  as  dis¬ 
solution  and  disturbance. 

Then,  if  D  =0,  the  system  consists  of  X,  C,  and  a  set  of 
ecological  response  functions  Re : 

X  «  Re  (C)  .  M 

If  D  ^  0,  we  must  consider  the  total  response  function  Rt* 

X  =  Rt  (C:D) .  (2) 

A  fundamental  problem  of  paleoclimatology  is  to  find  a  set  of 
■j. j-g j-jsf er  functions  0  such  that  C  can  be  estimated  given  X. 

C  =  0  (X)  .  <3> 

Three  points  are  worth  noting:  First,,  that  0  is  generally  obtained 
by  direct,  empirical  methods  and  not  by  inversion  of  Re  or  Rt* 
second,  that  C  is  a  matrix  generally  with  more  than  one  column; 
and  third,  that  when  D  ^  0  it  is  in  principle  still  possible  to 
derive  0  if  the  rank  of  X  is  greater  than  that  of  (C:D)  . 

The  three  elements  of  equation  (3)  are  defined  as  the  monitoring 
system.  The  system  includes  the  domain  X  and  the  range  C.  The  X  and 


C  used  to  derive  the  transfer  functions  are  the  calibration  data-set.  ? 
The  X  to  which  the  transfer  functions  are  applied  is  the  application 
data  set.  The  main  focus  of  the  conference  is  on  the  monitoring 
system  and  the  problems  of  obtaining  and  verifying  0.  Some  attention 
must,  however,  be  paid  to  the  problems  raised  by  equation  (2) ;  and 
the  opportunities,  such  as  those  explored  by  Wolfgang  Berger,  for 
finding  a  function  0^  to  estimate  D  given  X. 

i  '3 

i  j 

III.  OPTIMIZATION  PROBLEMS. 

In  order  to  derive  any  transfer  function,  a  number  of  choices 
must  be  made  among  alternative  procedures  (see  Appendix  1) .  The 
problem  therefore  arises  as  to  criteria  for  selecting  the  optimum 
set  of  procedures  for  a  given  monitoring  system.  We  wish  to  choose 
the  set  which  results  in  a  transfer  function  which  is  the  most 
robust  against  various  types  of  distortion;  the  most  precise  in 
terms  of  laboratory  error;  and  the  most  accurate.  Several  cate¬ 
gories  of  problems  can  be  recognized.  For  convenience  in  discussion, 
the  categories  are  given  letter  designations  and  the  problems  are 
given  a  single  numerical  code  and  a  brief  name. 

A.  Algebraic  problems .  Given  a  particular  domain  X  (in  biotic 
work,  a  set  of  taxonomic  categories  or  measures;  and  a  calibration 
set  of  samples) ,  the  following  problems  represent  choices  among 
alternate  procedures. 

. 

1.  Covariance-adjustment  problem.  Is  it  wise  to  "eliminate 

! 
i 


the  effect  of  latitude",  for  example. 
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2.  Data -sea ling  problems.  Various  data- transformation 
options,  including  those  spoken  of  as  “R-mode"  and  "Q-mode"  options, 
standardization,  row-normalizing,  and  others. 

3.  Statistical-technique  problem.  Generally  multivariate 
statistical  technique  or  series  of  techniques  should  be  used? 

a.  Eigenvector  problem.  Is  it  better  to  use  “raw  data" 
or  linear  combinations  defined  by  eigenfunctions? 

k.  .Rotation  problem.  If  eigenvectors  are  used,  should 
they  be  rotated?  Does  it  make  a  difference? 

c*  of  eigenvectors  (canonical  variates)  problem. 

How  many  of  the  calculated  eigenvectors  or  canonical  variates  derived 

from  X  or  (X:C)  should  be  retained  for  establishing  the  values  of  the 
transfer  functions? 

d.  Regression  problem.  What  procedure  is  best  for 
writing  0:  standard  least-squares  regression;  step-wise  regression 
with  elimination  of  less-useful  terms;  canonical  correlation;  etc.? 

e‘  Linearity  problem.  Should  0  be  assumed  linear  or 
non-linear?  What  is  the  nature  of  the  non-linearity?  If  the  non¬ 
linearity  consists  of  simple  squares  and  cross-products  of  X,  then 
the  problems  are  considered  under  this  heading.  If  a  more  complex 
form  of  non-linearity  is  involved,  the  problem  should  be  treated 
under  section  IV  below. 

4.  Confidence-ljgnit  problems .  How  are  these  to  be  assigned? 
Empirically  (i.e.,  with  an  independent  set  of  test  data);  or  theor¬ 
etically  (using  principles  of  statistical  inference  and  observations 
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on  the  calibration  data-set)?  What  good  are  they? 

B.  Domain  problems .  A  very  large  number  of  options  are  possible 
in  defining  the  domain  X,  i.e.,  in  specifying  the  kinds  of  variables 
to  measure  and  the  time-space  realm  in  which  to  measure  them.  These 
options  exist  for  both  biotic  and  isotopic  work. 

5.  Taxonomic  problems .  What  species?  what  morphotypes? 
what  combination  of  species?  what  combination  of  higher  taxa?  what 
habitat-selection?  etc. 

6.  Operational  problems .  What  census  size?  what  sieve 
fraction?  what  criterion  of  countable  specimen?  what  strew  procedure? 

7.  Realm  problem.  What  distribution  of  samples  in  time  and 
space  is  to  be  used  in  the  calibration  data-set?  This  is  a  very 
important  set  of  options,  with  implications  extending  into  questions 
posed  above  under  the  linearity  problem,  for  example,  and  several 
listed  below  as  fundamental  problems.  One  subproblem  is  to  evaluate 
the  method  used  by  Luz  (1973)  ,in  which  the  biological  measures  (the 
eigenvectors  derived  from  foraminifera  data  comprising  the  sea-bed 
calibration-data-set)  are  defined  by  a  data  set  of  coretops  and 
down-core  fossils  including  the  data  in  the  down-core  data-set  to  which 
the  transfer  functions  «ra  to  be  applied. 

8.  Inversion  problem.  Under  what  circumstances  does  0  not 
exist?  or  possess  limited  accuracy  (as  is  the  case  of  a  monospecific 
biota) ? 

9.  No-analogue  problem.  Occurs  when  fossil  values  of  cer¬ 
tain  taxa  exceed  the  modern  values  used  to  derive  the  transfer  func¬ 


tions  . 


c*  j*i!LDH.g.  problems .  The  selection  and  acquisition  of  valid 
estimates  of  the  physical  parameters  arrayed  in  the  matrix  C  are 
often  of  critical  importance. 

10.  Ground- truth  problem.  How  good  are  data  on  winter 
temperatures  for  the  Southern  Ocean,  for  example? 

Parameter-selection  problem.  Which  parameters  are 
ecologically  the  most  likely  to  influence  the  biotic  response,  or  to 
be  linearly-related  to  those  which  do?  Which  species  are  in  isotopic 
equilibrium?  Is  it  wise  to  use  such  parameters  as  atmospheric  pres- 
sure#  air-mass  frequency,  water  mass  frequency,  dynamic  height,  or 
sigma-t?  Should  we  use  standard  climatic  variables  or  use  linear 
combinations  or  complexes  of  these  variables? 

12  *  Parameter-independence  problem.  How  much  statistical 
independence  among  parameters  is  enough  to  provide  a  statistical 
basis  for  generating  transfer  functions  which  are  capable  of  giving 

independent  estimates  of,  say,  salinity  and  temperature;  or  rainfall 
and  temperature? 

F^i^er  problem.  Many  geological  samples  integrate  the 
climatic  inputs  over  many  centuries  or  more.  How  can- this  effect 
be  assessed? 

14.  Lag  problem.  Many  biological  and  sedimentary  monitor¬ 
ing  systems  exhibit  significant  lag  in  their  responses  to  climatic 
change.  How  can  this  effect  be  assessed?  How  do  we  deal  with  this 
problem  in  selecting  C  for  writing  0? 


15.  Seasonality  problem.  How  do  we  write  when  the 
calibration  data-set  for  a  single  hemisphere  crosses  the  thermal 
equator,  so  that  the  warmest  temperatures  occur  in  some  places  in 
August  and  in  others  in  February?  How  do  we  estimate,  say,  an 
August  sea-surface  temperature  map  for  the  Indian  Ocean,  18,000 
years  ago? 


IV.  FUNDAMENTAL  PROBLEMS. 

A  number  of  problems  require  consideration  at  a  different 
level  from  those  listed  under  section  III.  Some  represent  questions 
which  are  unanswerable  within  the  context  of  any  transfer  function 
based  on  a  single  biological  group  or  isotopic  system.  Others 
represent  a  questioning  of  assumptions  used  in  writing  transfer 
functions.  As  listed  below,  several  problems  overlap. 

16.  Validation  problem.  How  can  we  test  the  absolute 
accuracy  of  estimates  of  past  climate?  Four  approaches  have  been 
or  are  being  tried.  Any  others? 

a.  Direct  check .  In  tree-ring  work,  meteorological 
records  have  been  used  to  validate  estimates  of  the  recent  past. 
Where  else  can  this  be  done?  Do  these  tests  validate  estimates 


from  the  more  remote  past? 

b.  Comparison  of  two  or  more  independent  biologically- 
based  transfer  functions.  Radiolaria  and  forams,  pollen  and  tree- 


rings,  etc. 


estimates . 


d .  Concordant  estimates ,  i .  e . , 


pairs  with  overlapping 


confidence  intervals,  encourage  belief.  Like  discordant  isotopic 
ages,  discordant  estimates  (i.e.,  pairs  with  non-overlapping  con¬ 
fidence  intervals)  identify  important  methodological  problems  or 
geological  processes.  Discordancy  can  result  from  the  application 
of  two  transfer  functions  based  on  one  organism  group  cr  based  on 


more  than  one  group. 


.  consistency.  On  an  intuitive  basis ,  the 


areal  pattern  and  absolute  range  of  synoptic  maps  can  be  evaluated. 
The  advent  of  numerical  simulation  programs  offers  the  possibility 


of  making  such  intuitive  evaluations  both  more  rigorous,  and  more 


inclusive.  Ultimately,  proxy  data  from  tree-rings,  pollen,  plank 
ton,  isotopes,  and  glacial  geology  must  fit  dynamically. 

17.  Distortion  problem  (and  opportunity) .  How  effectively 
can  transfer  functions  eliminate  or  surpiress  the  effects  of  various 
processes  which  tend  to  distort  the  paleoclimatic  signal?  The 
other  side  of  this  coin  is  the  opportunity  to  obtain  useful  estimate 
of  the  intensity  of  the  distortion-producing  processes  themselves. 

Such  processes  include: 

a.  Dissolution  processes 

b.  Bottom-transport  processes. 

c.  Cultural  interference 

18 __  Quality— control  problem .  Can  &  priori  criteria  be  formu¬ 
lated  that  can  identify  samples  on  which  a  given  transfer  function 
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will  give  invalid  estimates?  For  instance,  the  criteria  can  be 
red  flags  of  the  form.  "If  species  j  exceeds  abundance-level  q, 

watch  out! " 

19.  Evolution  problem.  If  elements  of  the  biota  have  evolved 
since  the  fossil  deposit  was  formed,  the  calibration  and  application 
data-sets  are  non-homogeneous ,  and  trouble  lurks.  Analagous  pro¬ 
blems  occur  where  the  calibration  data-set  is  not  wide  enough  to 
include  the  present  range  of  fossil  species  (see. 9,  no-analogue 

problem) . 

20 .  Ecosystem-stability  problem.  What  happens  if  links  .be¬ 
tween  various  physical  elements  of  the  modern  ecosystem  (as  cap¬ 
tured  in  the  calibration  data-set)  become  uncoupled?  If  such  an 
uncoupling  occurs,  does  it  result  in  an  uncoupling  between  the 
biotic  and  physical  realms?  Do  species-responses  —  as  captured 
in  augmented  and  inverse  form  in  0  —  reflect  the  action  of  in¬ 
dividual  physical  parameters  themselves;  or  do  these  responses 
reflect  the  influence  of  these  parameters  as  they  are  locked  to¬ 
gether  in  a  particular  pattern  within  the  realm  of  an  air-  or 
water-mass? 

21.  Bio-thanatocoenosis  problem.  Are  increases  in  accuracy 
and  scope  awaiting  the  mathematical  modeling  of  isotopic  and  separate 
biological  steps  in  the  environment— plankton— fossil  and  environment— 

vegetation-pollen  chains? 
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22.  Transfer- function  structure  problem.  What  form  of 

0  should  be  assumed?  If  non-linear,  what  kind  of  non-linearity 
(see  3e)?  To  what  extent  is  the  form  dependent  on  data-scaling 
and  the  use  of  eigenvector  transformations  (2  and  3a)?  Do  dif¬ 
ferent  biological  groups  require  different  kinds  of  transfer 
functions? 

V.  TERMINOLOGICAL  PROBLEMS 

23.  Is  the  term  transfer  function  a  good  one.  What  about 
calibration  function? 

24.  Faunal  Index  vs.  faunal  temperature-estimate,  paleo- 
temperature  vs.  isotopic  temperature-estimate. 
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CORRESPONDING  PATTERNS  OF  CONTEMPORARY 
POLLEN  AND  VEGETATION  IN  CENTRAL  NORTH  AMERICA 

T.  WEBB  III,  Department  of  Geological  Sciences,  Brown  University, 
Providence,  R.  I.,  02912,  U.  S.  h. 
and 

J.  II.  McAndrews,  Department  of  Geology,  Royal  Ontario  Museum, 

Toronto,  Ontario,  Canada 

The  increased  use  of  modern  pollen  as  a  basis  for  interpreting 
fossil  data  requires  compilation  of  the  modern  data  in  readily 
accessible  form,  such  as  maps  of  the  relative  frequencies  of  indi¬ 
vidual  pollen  types.  This  paper  undertakes  this  task  and  presents 
maps  incorporating  over  600  pollen  samples  distributed  throughout 
central  North  America  (35-70°N,  75-110° W) .  Only  data  published 
after  1960  are  included,  and  data  from  50  sites  are  presented  for 
the  first  time.  The  maps  show  well-marked  changes  in  the  pollen 
from  plant  formation  to  plant  formation.  For  example,  peaks  in 
the  values  of  sedge  and  birch  occur  in  the  tundra;  high  values  of 
spruce  appear  in  the  boreal  forest;  high  values  of  pine  and  birch 
appear  in  the  mixed  conifer-hardwood  forests;  high  values  of  oak, 
elm,  and  hickory  occur  in  the  deciduous  forest;  and  high  values 
of  herb  pollen  appear  in  the  prairies.  Where  data  permit,  finer 
scale  patterns  in  the  vegetation  are  also  evident,  such  as  the 
difference  between  northern  hardwood  and  pine  forests.  Principal 
component  analysis  of  the  pollen  data  helps  to  illustrate  the  pat¬ 
terns  of  covarying  pollen  types  extant  in  the  data.  The  pollen 
maps  also  indicate  where  additional  samples  need  to  be  taken,  such 
as  north  of  65°N,  south  of  Hudson's  Bay,  and  in  the  plains  states. 
Joining  the  data  from  central  North  America  together  with  the  data 
from  eastern  North  America  compiled  by  R.  B.  Davis  and  Webb  pro¬ 
vides  the  initial  core  of  data  and  maps  needed  to  produce  an  "Atlas 
of  Pollen  Maps  from  Eastern  North  America." 
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-  — s:=  wttJssx  " 

“kcs  “ 

and 

T.  V'EBB  HI,  Dept,  of  Geological  Sciences,  Brown  Univ..  Prov.,  R.I.. 

U.S.A. 

,  _ -.j. r->  f  modern  do!  len  and  the 

An  understanding  of  the  relationship  e  ween  sp  a|ds  interpretation  of 

listribution  of  vegetation  on  a  c°"+'nen+a'  sS“',S, vidua  I  genera  were  constructed 
I a+e  Quaternary  assemblages  of  po  •  P  55  prev-,0usly  unpublished  samples,  eas 
from  400  samples  of  modern  pollen,  lncludin9J,  Isonoils  on  these  maps  clearly 

o™0°“n  North  America.  The  Iso  requency  odours  o^sopo^  ^  da+ 

reflect  vegetational  patterns  both  on  the  Particular  I  y  clear  for  several  of  the 

quite!  on  a  finer  scale  Correspondence^^  par^ar^  and  C^) 

+nxa  (e.q.  Cyperaceae,  P i cea_,  IP i nu_s_,  ,  ..p  assemblages  of  fossil  pollen, 

frequently  used  for  climatic  interpre  a  bers  Qf  surface  samples  in  areas  poorly 

The'maps  'also  Indicate  need  for  ncreased  numbers  ‘ yfj sort  re.a- 

sampled  at  present.  A  more  complete  und. ef+auant?ties  of  each  important  taxon  in 
ttonshi p  awaits  the  aval M -bill* foi  ^ of  fossil  pollen  however, 
+h/=.  v/Afip+at i on •  Potential  smbigul  y  . »  recorded  on  a  relative 

pot len^concentratto^decreos^nghnorthward  from  the  boreal  forest  into  the  region 
of  tundra. 
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ABSTRACT 

Emerged  coral-reef  terraces  on  the  Huon  Peninsula  in  Hew 
Guinea  v/ere  reported  in  e  reconnaissance  dating  study  by  Veeh  and 
Chappell  (1970).  Age  definition  achieved  was  not  good  for  several 
important  terraces,  and  we  report  here  a  series  of  new  Th2!VU254  dates, 
which  further  clarify  the  history  of  late  Quaternary  eustatic  sea-level 
fluctuations.  More  than  20  reef  conplexes  are  present,  ranging  well 
beyond  250,000  years  old  :  we  are  concerned  with  the  seven  lowest 
complexes.  Major  reef-building  episodes  dated  by  Th2 3 °/U* ***  are  reef 
complex  I  at  5  to  9  ka  (ka  =  1000  yrs),  r.c.  Illb  at  41  ka  (4  dates) 
r.c.  IV  at  60  ka  (4  dates),  r.c.  V  at  30  ka  (2  dates),  r.c.  VI  at 
106  ka  (2  dates),  and  r.c.  VII  at  118  to  142  ka.  Complex  II  was 
previously  dated  by  C14  at  28  ka  :  this  age  has  not  yet  been  confirmed. 
The  reef  crests  were  built  during  or  immediately  before  intervals  of 
sea-level  maxima,  when  rates  of  rising  sea  level  and  tectonic  uplift 
briefly  coincided.  The  culmination  of  each  reef-building  episode  was 
only  a  few  thousand  years  in  duration,  and  multiple  dates  from  the 

same  reef  complex  generally  group  within  the  statistical  errors  of 

* 

the  individual  dates. 

Several  methods  can  to  us.-d  to  estimate  the  altitude  of  each 
sea-level  maximum  relative  ;o  present  sea  level.  The  least  complicated 
is  to  calculate  m:*an  tectonic  uplift  rate  for  each  profile  of  the 
terraces,  and  use  the  mean  rate  to  calculate  the  tectonic  displacement 
of  each  dated  reef  complex  on  that  profile.  The  difference  between 
the  present  altitude  of  a  reef  complex  and  its  calculated  tectonic 
uplift  gives  the  palcosea  level  at  the  time  the  reef  grew.  He 
estimate  uplift  rates  for  six  surveyed  sections  by  calibrating  against 
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published  paleosea  level  estimates  from  Barbados  and  elsewhere,  viz. * 

125  ka,  paleosea  at  +6  m;  103  ka,  -  15  m;  82  ka,  -  13  m.  For 
each  section  the  individual  uplift  rates  for  reefs  V,  VI,  and  V 1 1 b 
are  within  5 %  of  their  section  means.  Using  the  mean  rates,  paleosea 
level  estimates  for  reef  crests  II,  IIIB  and  IV  are  made  for  each 
section.  Consistency  of  estimates  between  sections  is  good,  giving 
-28  m  for  the  60  ka  paleosea  level ,  around  -38  m  for  the  42  ka 
level,  and  -41  m  for  the  28  ka  level  (age  unconfirmed).  Using  the 
mean  uplift  rates,  the  82  ka  and  103  ka  paleosea  levels  are  also  estimated 
for  each  section  :  all  individual  estimates  are  plotted  graphically, 
and  a  sea  level  curve  drawn.  The  reef  stratigraphy  indicates  sea- 
level  lowerings  between  each  dated  reef  crest  :  the  crests  probably 
represent  the  interstadials  of  the  Wisconsin  (Wlirm,  VJeichsel) 

Glaciation,  and  intervening  lower  levels  correpond  to  stadials.  Since 
the  last  time  of  eustatic  sea  level  higher  than  the  present  (about 
125  ka),  5  sea-level  maxima  occurred  at  roughly  20  ka  intervals,  none 
being  as  high  as  the  present. 


Reef  Complex  I 


?S  i  -i 


111-era 
Gor.ias  fcraa  reiiformis 
lnoobora  miaroconcj 


Loo tor j 


i?.  stallieei-a 


3.13+  .Oo 

1.32+  .02 

0.17+  .01 

21-^-Z 

2.63  +  .Op 

1.12  +  ,C? 

0.037+  .005 

9  +  0.5 

2.48  +  .05 

1.15  +  -'ll 

0.051  +  .00'; 

5  +  0.4 

2.64  +  .05 

1.12  +  .02 

0.037  +  .005 

9  +0.5 

3.iq  +  .o: 

1.11  +  .02 

0  035  +  .00’+ 

9  +  0.4 

Pc.via  sm. 

Reef  Complex  II 
2.59  +  .08 

1 .11  +  .03 

0.031  +  .002 

T.erjtoria  phry^ia 

2.80  +  .08 

1.12  +  .03 

0.05’+  +  .003 

Reef 

Complex  III 

.  • 

Pavia  sneciasa 

3.01  +  .06  1.13  +  .02 

0;23 .+  .02  ; 

'35  +  3 

ICrir-cahora  exes  a 

2.70  *  .05  1.13  +  .02 

0:32  +  .02 

41  +3 

Genies trea  parvistella 

2.54  +  .05  1.09  +  .02 

0.32  +  .02 

42+3 

Iraboohvllia  coryr.bost 

3.58  +  .07  1.12  +  .02 

O.O69+  *004 

7  +  0.4 

Symn'ayllia  r.obilis 

3.23  +  .10  1.10  +  .03 

0.32  +  .02 

41  +  3 

- 

Pavia  stellieera 

Reef  Complex  IV 

2.39  +  .05 

1.11  +  .02 

0.1*2  +  .02 

58+4 

Fa  via  pallida 

2.37  +  .05 

1.11  +  .02 

0.43  +  .02 

61  +  4 

Acropora  sp. 

4. Co  +  .08 

•  1.15  +  .02 

0.33  +  .02 

45  +  3 

Favia  pallida 

3.12  +  .06 

1.11  +  .02 

0.4l  +  .02 

57  +  4 

Hydnobhora  microcodes  2.92  +  .06 

1.13  +  .02 

0.46  +  .02 

66+4 

Reef  Complex  V 

■ 

Haty^yra  lamellir-a 

2.63  +  .05 

1.13  +  .02 

0.43  +  .02 

60  +4. 

Porites  lutea 

3.49  +  *07 

1.10  +  .02 

0.54  +  .01 

84+2 

Ooniastrea  oectinata 

2.9'+  .06 

1.11  +  .02 

0.55  +  .02 

85+4 

Reef  Complex  vj. 

Pavia  soeciosa  2.73  +  .Co 

1.12  +  .02 

0.63  +  .03 

106+7 

Kvir.odr.ora  r.lerocoros  3.45  +  -Oi 

1.09  +  .02 

0.63  +  .02 

107  +  6 

Porito;  lutaa 


2  9'/  +  .Co  ?.3B  4  .08  0.74  4  .02  14?.  +  8 
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TABLE  2.  MINERALOGY  AND  MINOR-ELEMENT  CHEMISTRY  OF  LATE  QUATERNARY  CORAL  SAMPLES 


FROM  THE  HUON  PENINSULA,  NEW  GUINEA  (SEE  TEXT  FOR  DISCUSSION). 
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Table  2  (Contd.) 


39a 

2,3,c 

.90 

- 

.10 

6.20 

1.50 

39b 

1 

.99 

- 

.01 

6.70 

1.36 

39b 

1 

.99 

- 

trace 

7.10 

1.07 

1 

1.00 

- 

- 

6.40 

1.24 

41 

2, a  ,b 

.96 

- 

.04 

6.20 

1.22 

42 

1 

.99 

- 

trace 

7.60 

2.06 

43 

2 

.92 

.04 

.04 

7.60 

2.40 

44 

2,3,a,c 

.94 

- 

.06 

6.44 

1.92 

45 

1 

1.0 

- 

- 

7.56 

1.29 

46 

2,a,b 

.95 

- 

.05 

6.70 

!  2.04 

1  dated  (see  table  1  for  additional  data) 

2  rejected  for  dating 

3  thin  section  prepared 

a  void-filling  crystals 

b  void-filling  micrite 

c  calcite  crust  in  corallites 

d  micrite  coating  in  corallites 

e  borings 

f  hairline  of  clear  calcite  crust 
g  internal  sediment  and  aragonite  needles 
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TABLE  3.  MEASURE  OF  REEF  CREST  ELEVATIONS  (m) ,  AND 
UPLIFT  RATES  (m/103  YRS) ,  CALCULATED  USING 
ASSUMPTIONS  GIVEN  IN  TEXT: 


TERRACE 

Age, 

Initial 

Elevation 

KANZARUA 

BLUCHER 

KWAMBU 

KAMA 

SAMBERO 

KAMBIM 

VI  lb 

124  ka, 

+6  m 

330, 

2.62 

280, 

2.21 

215, 

1.69 

160, 

1.24 

150,' 

1.16 

120,'. 

•  *91 

VI 

103  ka, 

-15  m 

250, 

2.57 

215, 

2.23 

160, 

1.69 

115, 

1.26 

no, 

1.21 

93 

1.05 

V 

82  ka, 

-13  m 

190, 

2.48 

155, 

2.05 

117, 

1.58 

90, 

1.25 

80, 

1.13 

60 

.89 

mean  rate 

2.56 

2.16 

1.65 

1.25 

1.17 

.94 

predicted  height  of 
reef  I  (m)  * 

10.2 

8.6 

6.6 

5.0 

4.7 

3.8 

actual  height  of 
reef  I 

15 

10 

6 

5 

5 

2.5 

*  calculated  as 

uplift 

rate  x  4 

ka;  see 

text. 

TERRACE 
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TABLE  4.  MEASURE  OF  REEF  CREST  ELEVATIONS  (m) ,  AND 
RECONSTRUCTED  PALEOSEA  LEVELS  (m),  USING 
MEAN  UPLIFT  RATES  FROM  TABLE  3,  FOR 
TERRACES  II  TO  IV. 


Age, 

x  1000  yr 


< 

ZD 

CZ 

CsZ 

UJ 

*=c 

zsz 

N4 

CJ 

zz 

ZD 

<£ 

-_l 

CO 

ZD 

CO 

s: 


o 

C£ 

s 

Ll) 

i — * 

CQ 

s: 

CO 

s: 

<c 

<c 

< 

z 

to 

1Z 

Mean 

(m) 


IV  60  125, 

-29 


70,  48, 

-29  -27 


28,. 

-28  ■  ' -28 


Ilia  * 

50  max. 

40  min. 

90 

-38 

-12 

65 

-43 

-21 

42 

-42 

-26 

"  5 

*  » 

see 

text 

Illb 

40 

70, 

-32 

41, 

-45 

28, 

-38 

10, 

-40 

10, 

-37 

"  y 

-38 

II 

28 

30, 

-42 

18, 

-42 

7, 

-40 

"  1 

“  1 

"  i 

-41 

NOTES  :  *  As  Ilia  has  not  been  dated  directly,  paleosea  levels 

are  calculated  for  upper  and  lower  age  limits  -  see 
text. 

-,  Terrace  crest  uncertain;  or  no  satisfactory  height 
measurement. 
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CAPTIONS  FOR  FIGURES  AND  TABLES 


Figure  1.  Location  and  generalized  physiography  of  the  Muon  Peninsula 
New  Guinea.  Spot  heights  in  meters;  rainfall  in  inches. 

Figure  2.  Three  aerial  views  of  the  Huon  Peninsula  emerged  reef  terraces 
,  (Chappell,  1972).  Variations  in  elevations  of  reef  complex 

VjjH —  III  and  VII  are  a  guide  to  the  differential  uplift  of  the 

region.  Note  :  altitudes  on  uppermost  photo  are  barometric 
estimates  only;  other  heights  are  surveyed. 

Figure  3.  Models  of  the  facies  geometries  in  barrier  and  fringing  reefs 
showing  the  changes  associated  with  submergence  and  emergence* 
(a)  Landward-dipping  internal  facies  boundaries  associated  with 
submergence,  (b)  Changes  in  reef  morphology  during  submerqence 
from  (1)  fringing  reef,  (2),  (31  barrier  reef  with  lagoon  and 
fringing  lagoon-shore  reef,  (1)  barrier  reef  with  lagoon  and 
patch  reefs.  (c)  Cliffs,  notches,  and  fringing  reef 
associated  with  emergence  followed  by  relative  stillstand. 

Figure  4.  Six  surveyed  profiles  of  the  Huon  Peninsula  emerged  reef  • 
terraces.  Reef  complexes  I  to  VII  are  shown  by  Roman 

numerals.  Sample  numbers  are  shown  as  near  as  possible  to 

their  collection  localities. 

Figure  5.  Late  Quaternary  paleosealevels  based  on  estimates  from  New 
Guinea  and  elsewhere. 


Table  1.  Sample  numbers,  taxonomy,  radiochemistry,  and  ages  of  late 

Quaternary  coral  samples  from  the  Huon  Peninsula,  New  Guinea. 

Table  2  Mineralogy  and 'minor-el ement  chemistry  of  late  Quaternary 
coral  samples  from  the  Huon  Peninsula,  New  Guinea. 

’  Table  3.  Measured  elevations  and  calculated  uplift  rates  for  reefs 
V,  VI,  VI lb,  for  six  surveyed  traverses. 

Table  4.  Measured  elevations  and  paleosea  level  estimates,  using 

uplift  rates  from  Table  3,  for  reefs  II,  Illb,  IHa,  IV  for 
six  traverses.  * 


VI  Reel  complex  number  j , 

*  Theodolite  survey  point  ^  1  * 

32  Field  sample  from  traverse,  as  shown, 

4l'  Field  sample  from  near,  but  not  on,  traverse 


Example 

Time  interval  Rate  of  Uplift 


ADC 


T1 

to 

T2 

2.0 

1.0 

0.5 

T2 

to 

T3 

1.0 

0.5 

0.25 

T3 

to 

T4 

4.0 

2.0 

1.0 

R  -  2.33 

1.167 

0.583 

Input  to  model 

Estimated  from  rate  differences  I 

R.  =  c  •  5  ,  and  likewise 

A  A 

n 

for  traverses  B  and  C 

-  v„-  -  V 

likewise  for  (R^  -  Rc 

and 

ave. 

)  and  (Rp  -  Rc) 

(RA  -  Vn 

(ra  -  V 

<»B  -  »c' 

Time  interval  c 

c 

c 

c 

T  to  T2  .857 

.857 

.857 

.857 

T2  to  T3  .429 

.429 

.429 

.429 

T,  to  T  ,  1.715 

3  4 

1.715 

1.715 

1.715 

c  =  1.000 

1.000 

1.000 

1.000 

Thus,  a  set  of  coefficients  (c)  constituting  a  variable  uplift  model 

can  be  estimated  from  the  rate  differences  (R^  -  R  )  ...  Knowing  R  .. 

A  B  n  A 

from  rhe  age  and  elevation  of  any  one  terrace  of  known  sea  level, 
estimates  of  all  other  sea  levels  (n)  can  be  calculated  by  application 
of  the  uplift  model  to'  each  traverse. 

R.  K.  Matthews,  March  1974 
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R«produc«dlrom 

b«it  «v»il«bU  copy* 

SEA  LEVEL  DYNAMICS  AT  THE  END  OF  TIIE  LAST  INTERGLACIAL: 
INFEEBsJCiSS  FROM  125,000  YEAH  01.12  REEF,  BARBADOS,  W.  I. 

Many ion,  Stephen,  and  Matthews,  R.  K. ,  CLIMAP,  Dept. 

Cool.  Sci.,  Brown  University,  Providence,  R.  I.  02912 

Spatial,  relationships  of  coral  tuxa  within  the  125,000  year 
old  Barbados  III  reef  indicate  sea  level  during  the  Pan¬ 
gs!  on  maintained  its  maximum  position  for  less  than  5,000 
years,  before  abruptly  and  rapidly  receding.  Development  of 
a  reef  facies  geometry  was  controlled  by  the  rate  of  (1) 
tectonic  uplift,  (2)  reef  carbonate  production,  and  (3)  sea 
level  fluctuation.  Analyses  of  cross-sections  and  borehole 
samples  of  the  reef  have  been  used  to  interpret  the  follow¬ 
ing  sea  level  history.  Early  phases  of  reef  accretion  are 
marked  by  nearly  vertical  contacts  between  adjacent  facies 
indicating  a  rapid  sea  level  rise.  Additional  diagenetic 
evidence  from  borehole  samples  indicates  this  rise  started 
at  least  20  m  below  present  sea  level.  Later  phnses  are 
marked  by  a  maximum  of  40  m  of  reef  progradation  with  sen 
level  at  or  near  its  highest  elevation  for  no  more  than 
5,000  years.  Finally  the  reef  crest  facies  moved  down- 
,  slope  as  reef  construction  was  progressively  halted  by  a 
sea  level  lowering  of  5-10  m/1,000  yrs.  The  rcite  of  sea 
level  lowering  indicated  by  this  study  implies  that  post 
Panyaion  glacial  growth  was  rapid. 
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Chnirn&n: 


J/VMES  D.  HAYS  (Lnr.ont-Doherty  Go ological 
Observatory,  Palisades,  Kdvj  York)  end 
J.  I1RRAY  MITCKEU,  JR.  (s.Utt.  Environmental 
Data  Service,  Sil"cr  Spring,  Maryland) 


^Kurr^roa:  se*£  sbsST 

clal  Ira  durlngEtHetlastCnilllonkyca«\ad  th*3'  "d"6  naJ°r  advances  ond  retrcats  of  gla- 
nagnltude  This  13  (n  'l  l  ?"d  th?5e  “dvancc3  znS  retreats  were  of  similar 

cate  slmtarlcl  Snk  ET2m.~  5  ,*  V ^  *laclal  deposlcs  on  ^nd  •*»««•  iHi- 

the  slaclal^Lrglac^^c^  freJear  r«f  o«S*d  1°'  th°  laSt  5°°*000  >eats* 

oral  simultaneous  ilm  of  "  “KV  $  afd^S  £££*  ‘’nV1^  *“ 

nff  -s  noavccrte  (the  ^  30'0°0  -;)°thes;es  0° h0T  ;e  It  It YZ*Z£2Z 

ftnde  occaas  in  the  N.  Hemisphere  lag  the  thermal  response  in  the  tropics.  2 
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roi^rlTe'^£r?-' ^^U-r-1-^-^^s-^501O00J(ears.  Pleistocene  sea  level  is  largely 
ro....ol.ed  by  continental  ice  volume.  Coral  reef  sequences  and  intervening  subacrial  * 

luH  riLrifi?'  I  hfS,C°Iy  05  ?Sa‘1<!Vel  «»«tuatUn»  that  can  be  dated  by  Th-230 
jugnout  the  last  full  glacial-interglacial  cycle.  Interglacial  sea  lev-1  125  nnn 

reached  n8"nQ  belowCSCnt  1  ’  ?*“?  hct',°en  tha  125.000  and  105.000  high  stands  ‘ 

I^us  ‘the  d«“ioplv  fha  ^  1  f\Vl,thl?  5*03°  »“**  °f  thc  p«v*°-  hl«h  «a"d- 
quite  lapidlyf  ?  ‘  C°ntlnental  ^^“Ts  can  buiid  tn  substantial  proportions 
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R.  K,  Matthews 
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ni  Tgr:ga;r  ~C~ST5-C-— --  ~h-lr-^-bal  Distribution  18,000  yrarc  Ago  and  Subsequent 

Antarctica  by  6.000-10  000  years  ago  ’  A  Darti0"!  S*0’  °f  111  *Ce  shcct3  outsitle  of 
rapid  disintegration  of  .hJT ?  *  particularly  spectacular  event  involved  the 

8.000  yearl  Z?  mnfr  oclio^' Tf  J  ^  P°“lon  of  the  J-n““ntlde  Ice  Sheet  about 

occurred  at  2,500  year  Intervals  The'Ent'Am’"1?!'  *upe^DP03ed  on  the3e  naJ°r  events, 
cntal  base.  h«  reLl^ed  rll^ly  «able.  bu  ^e  ^Ct%WM- -h  r“CS  3  C°ntln* 

SSaSjfesSSr 
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Ji^Spa££^tjern_oLCli^te  Change  in  the  Atlantic  Ocean.  Per  sediments  beyond  the 
f"n  varla'tionrin°L  COrtelatdva  h°.rlzona  in  deep-sea  sediments  are  interpreted 

wa-ks  llie  1  a  ►  *  fl0ra'  aad  scdlncnt  c'l^aistry.  Thc  resulting  correlation  nct- 
U.ks  have  been  used  to  construct  a  transect  of  the  Atlantic,  which  shows  vateroass- 
nlgrations  during  the  last  200.000  years  across  rore  than  20*  of  latitude.  At  peak  gla- 

“ter  ?°T  ‘TV"  42’JI  U,I"°  a"  °!>rupt,  east-west  frontal  system  sep¬ 
arate  an  expanded  cyclonic  subpolar  gyre  from  the  subtropical  gyre.  Comprehensive  Eto^ 
g  a?.'  c  core  coverage  indicates  that  during  the  last  degiaclation  (18,000-6,000  years  ago) 

'  n,  c  a  9°Ioor:  Ir’old  ItS  ‘=lacU2  ail£nment.  in  the’sfbp^"  H.£ ^ 

ocenU  floni  alo^  PCrDltS  3  s5m0?tla  reconstruction  which  defines  an 

and  2o!i>IlIeL  flum  rh  ^  a?Pa“a'ly  «aatud  by  a  maximum  wind-stress  axis  anchored  along 
^  fru”  the  southcrn  “rKin  of  the  Laurentide  ice  sheet.  Similarly  the 
r  a  r11  0f  thc  Ice-^c  P»lar  front  (18.000  years  ago)  sfg- 

LZA  I  ^  U-1S  ln  "t3bll3hcd  b>'  ^c  focusing  of  an  axis  of  rLimum  «st"ly 

ind.  ^.long  the  narrow  latitudinal  band  lying  between  the  southern  margins  of  the  ^ 
Laurent idc  .ind  Scsndiniivijn  lc»!  sheets* 
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)J.*i  1 .1  Oy«-nt  \ *.< »t h  e  r  j; , p JDu  cj  n?  r :i  _l c r  a - \  Ctj COn  Vi;. ir^_  /,  «-,o .  One  aiu  of  tho  C.t.lMAi* 

«  project  i‘i  to  recount net  Lite  surface  n:o iiogr.ijihy  of  tlm  world  ocean  for  particular 
tiru-.:  during  the  pa". t  700,030  yc.irs,  Planktonic  fo.isils  from  chtouor.ircft  I  graphic  samples 
iiio  processed  using  null  ivnr  lie  a  tv.mif  er~  fuactio.i  technique::  (calibrated  on  the  r.odcrn 
t;c.i  to  yield  ciiCiaUrA  of  seasonal  surface-water  tenpcr.i tur»*s  and  salinity.  The 

*  first  result  ir.  presented  here.  Considered  ns  an  atto.-aly  pattern  ne.tsm’cd  f roa  today's 
base-line,  tli**  therrnl  rospoase  of  the  sea-surf  act  16,000  yearns  ago  in  roughly  th“  j;.nr.?c 
l  ^  il.  U.  1  iA„.ill  vidv.  And  F.4ct.  lu  a*  a  t  ip  of  today's  seasonal  oin.,0.  Change?  m 

r.orv'  ti  id-la  tit  Usle  sites  reach  and  scv.ioc  lie;;  exceed  10°C.  Except  foe  up-welling  areas,  low 
latitude  change-?  are  on  the  order  of  2  C,  and  in  some  subtropical  areas  chow  no  change. 
Ocher  feature-,  are:  1|  th*  raided  0-4.1  *3 roreini  ihT  T*  Jrfs  «■  *****  l* 

Atlantic  and  Southnvn  Ocean,  but  not  iu  the  K.  Pacific;  2)  the  areal  expansion  of  sub¬ 
polar  water.;  in  all  oceans;  3)  the  stable  geciT.otry  of  t lie  subtropical  gyres  in  all  oceans; 
f*)  the  steepened  thermal  gradients  across  polar  fronts,  apparently  parking  the  axis  of 
ice-age  westerlies  in  both  hemispheres;  5)  aa  increase  iu  the  extent  of  sea  ice:  at  all 
star.c.i?.  i».  tie.  ti.  Atlantic,  and  especially  during  the  austral  summer  in  tne  Southern 
Ocean. 


Klements  of  0rd«»c  ir.  Ice-age  History: _ A  Markovian  View.  Evidence  of  order  is  found 

both  in  the  chronology  of  clioatlc  change  and  in  the  spatial  structure  of  ice-age  climate. 
A  quasi- periodic  glacial-interglacial  cycle  with  n  time  constant  of  about  100,000  years 
is  clearly  stamped  on  records  of  the  past  500,000  years.  During  the  peak  of  the  last 
ice-age  the  ocean  circulated  in  expanded  cyclonic  and  restricted  anticyclonic  gyres.  One 
hypothesis  is  that  the  first-order  response  of  the  system  is  partly  determined  by  the  val¬ 
ue  of  changing  boundary  conditions  (the  geographic  pattern  of  incoming  annual  and  season¬ 
al  radiation  as  fixed  by  the  earth's  orbital  parameters)  and  partly  by  the  state  of  the 
system  Itself  (including  the  global  ice-inventory).  Major  glaciations  arc  associated 
with  long  intervals  of  high  cerJdional-gradients  of  annual  incoming  radiation,  amplified 
by  seasonal  effects  on  the  N.  Kenispkare  snow  line.  They  are  conditioned  globally  by  geo- 
strophic  expansion  and  contraction  of  oceanic  gyres  which,  especially  in  the  N.  Atlantic, 
alter  the  pattern  of  poleward  energy-transport .  Superposed  on  this  first-orJer  pattern 
are  small  amplitude  fluctuations  of  unknown  origin,  such  as  those  of  a  minor  cycle  of 
about  2,500  years;  and  random  fluctuations.  In  this  view,  climate  is  a  non-s tat ionary 
Markoy  p UwUtq wUwH  1*  c  *r>*  t  tb  i  t  -r  1.  Ltk  *..d 

elements. 


Mump  deal  Simulation  of  Ice  uejCli_natje.  Using  the  paleoclimat ic  boundary  conditions 
assembled  by  CEIMAP,  on  experiment  is  in  progress  which  aims  to  simulate  numerically  the 
%  global  climate  which  accompanied  the  last  ice  age  some  18,000  years  ago.  The  first  re¬ 
construction  will  be  for  the  Northern  Hemisphere  summer.  The  distributions  of  sca-cur- 
fnced  temperature,  ice  extent  and  elevation,  and  surface  albedo,  together  with  a  sca- 
level  lowering  of  85m,  will  be  used  with  an  improved  version  of  the  two-level  atmospheric 
general  circulation  model.  An  integration  of  several  month's  length  is  envisaged  in 
order  to  characterise  the  mean  August  values  of  the  principal  climatic  elements,  such  as 
pressure,  temperature,  and  rainfall.  Verification  is  to  be  ruidu  against  independent  data 
on  continental  climate,  chiefly  pollen  records. 


Perr..^ctiv^s  on  Climatic  Change.  As  clearly  apparent  from  earlier  presentations  in 
this  session,  paleocliratology  is  well  on  its  way  at  last  to  achieving  the  status  of  a 
quantitative  science.  On  the  one  hand,  reconstructions  of  the  chronology  of  Quaternary 
climate  not  o,iy  are  capable  of  being  c/prsssed  in  phvslctil  terns  readily  interpre table 
by  tho  atmospheric  scientist  but  enjoy  absclute  dating  controls  sufficient  to  permit 
reliable  correlations  with  datable  cnvironxental  events,  thereby  facilitating  the 
intelligent  investigation  of  causal  mechanisms  of  the  ice  ages.  On  the  other  hand,  the 
realistic  nur.orlcal  simulation  of  past  clinates  by  means  of  general  circulation  modeling 
experiments  is  fast  becoming  a  realizable  goal,  with  ran/  encouraging  results  already 
achieved  and  little  reason  to  doubt  that  both  climate  theory*  and  computer  technology  will 
suffice  in  the  rather  near  future  to  pin  down  both  the  rechanlsn  and  the  root  cause  (or 
causes)  of  the  grand  glacial-interglacial  succession  of  the  Quaternary.  Against  this 
auspicious  background,  the  present  state  of  our  understanding  of  the  ice-age  problem, 
and  that  of  the  related  problem  of  long-range  climate  prediction,  will  be  briefly 
assessed. 
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